Correlation functions after a quantum quench in a spin-chain in the pre-thermahzed 

regime 
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Results are presented for a two-point correlation function of a spin-chain after a quantum quench 
for an intermediate time regime where inelastic effects are weak. A Callan-Symanzik like equation 
for the correlation function is explicitly constructed which is used to show the appearance of two 
scaling regimes. One is for spatial separations within a light-cone, and the second is for spatial 
separations on the light cone. In both these regimes, the correlation function is found to decay with 
power-laws with new nonequilibrium exponents that differ from those in equilibrium, as well as from 
those obtained from a quadratic Luttinger liquid theory. 
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Motivated by experiments in cold atomic gases [l| 
and ultra-fast spectroscopy of strongly correlated ma- 
terials 0, the nonequilibrium dynamics of interact- 
ing quantum systems has become a topic at the fore- 
front of research. In this context, dynamics arising 
due to a quantum quench, where a system is prepared 
in the ground state of an initial Hamiltonian Hi, and 
then time-evolved with respect to a final Hamiltonian 
Hf, is of particular interest because of its potential for 
addressing several fundamental questions [4]. Some of 
these are, the mechanisms and time-scales for thermal- 
ization (Hl-liol , the possibility of an intermediate time pre- 
thermalized regime [lll-fis}. dynamical phase transitions 
associated with non-analytic behavior during the time 
evolution dynamics of integrable models flB- 

2l| and the possibility of describing their steady-state in 
terms of a generalized Gibbs ensemble (GGE) Yet 
another important question, which is related to the topic 
of this paper, is the a ppe arance of universal behavior 
in the dynamics [1, [23l. l24l| , with the possibility of cap- 
turing such a behavior in a renormalization group (RG) 
approach, even though the system is out of equilibrium. 

We study a quantum quench in a 
generic one-dimensional spin-chain H = 

-\- n.n.r 



where 



n.n.n. denote additional next-nearest-neighbor cou- 
plings. We study the dynamics in a continuum theory 
by taking into account the effect of the leading irrelevant 
operator. In particular we study the time-evolution of 
the two-point correlation function of the staggered spin 
component R{j) = (— 1)^ (S'qS'|) by explicitly deriving 
a Callan-Symanzik (CS) like equation for R, thus 
revealing conditions under which scaling holds in the 
nonequilibrium problem. We show that near the critical 
point, where the interaction is marginally irrelevant, 
logarithmic corrections arise (as in equilibrium 25- 27|). 
however the nature of the logarithmic corrections are 
qualitatively different for the nonequilibrium problem. 

We study a quench where initially the system is in the 



ground-state of a Luttinger liquid [2? 



K„ {7TU{x)f + (9,0(.T))2 
Ao 



(1) 



—dx4'/Tr represents the density, 11 is the variable canoni- 
cally conjugate to (j), Kg is the dimensionless interaction 
parameter, and uq is the velocity of the sound modes. 
The system is driven out of equilibrium via an interac- 
tion quench at t=0 from Kq — > K, with the leading ir- 
relevant operator corresponding to a commensurate lat- 
tice Vsg also switched on suddenly, at the same time as 
the quench. This triggers non-trivial time-evolution from 
t > due to Hf = HfQ + Vsg, where Hfo = Hi{Ko 
K,uo — >• u), Vsg — — / da:cos(70), and A — u/a is a 
short-distance cut-off. We assume that the quench pre- 
serves Galilean invariance so that uK=uoKq. We study 
the time-evolution of the equal time two-point correlation 
function of the staggered spin component. 



R{xit, X2t) — 4(cos 



J(j){Xit)\ ("f(j){x2t) 



cos 



(2) 



2 J — y 2 

Denoting 0{xt) = 2 cos ^^^^^^^^ j , R may be written as 
a Keldysh path-integral representing the time-evolution 
from the initial pure state (hence an initial density 
matrix p = \il>i) corresponding to the ground state of 
Hi, R{xit,X2t) = Tr[p{t)O{xi0)O{x20)] which may be 

written as R = Tr e-^"i^\'ipi){tpi\e'"f^d{xi)d{x2) = 

^V[(t>cu4>q]e'-''^°+^'^^di{xit)di{x2t), where 6/ is the 
operator in the interaction representation of -ff/o, S|n de- 
scribes the nonequilibrium Luttinger liquid (.9=0) [29|, 

So^t; I dxi I dX2 I dti I dt2{(l)cl{l) (t>q{^)) 



G-^\l,2) -[G^IG^G;^!] (1,2) 



0,(2) 



(3) 



a.nd Ssg^^JZ^dx.Jl^dh [cos(70_(l)) - cos(7</'+(I))]. 
Above l(2)=(a;i(2),ti(2)), with -/+ rep- 

resenting fields that are time/anti-time ordered on the 
Keldysh contour (soj . 
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In a global quench like the one we study, the system is 
translationally invariant in space, so that R{xit,X2t) = 
R{xi — X2,t). We define the following three exponents 
which play an important role in the dynamics. 



K — - ■ K — — 



Ko 1 



Ktr = -^K. 1 



^0 



(4) 



K(,q governs the power-law decay in the ground state of 
i?yo, Kneq determines the power-law decay at long times 
after an interaction quench in the Luttinger liquid 8, 31], 
Ktr determines the crossover from a short-time behavior 
determined primarily by the initial Luttinger parameter 
Kneq + Ktr — / ° , to the loug time behavior deter- 
mined by Kneq iM- The role played by these exponents 



is apparent in R which for 5 = and in the scaling limit 
(setting u—1 and r — \xi — X2\, t are in units of A) r,t 3> 1 
and far from the light cone (|r ± t| :$> 1) is |15l l31j. 

1/4 



i?(°)(r,i) = 



1 



(2t)2 



-Ktr 



(5) 



R shows the "horizon-effect" where outside the light-cone 
r ^ 2t, the correlator depends only on the initial wave- 
function (and hence the initial Luttinger parameter Kq) 
albeit with a time-dependent prefactor t^*' /^. However 
within the light-cone r <C 2t, the correlator reaches a 
steady-state characterized by the nonequlibrium expo- 
nent Krten- Morcovcr, on the light-cone R^^\r = 2t) ~ 

1/4 



^neq- 



-Ktr 



It is interesting to contrast this 



behavior with that of a correlator in a quench from an 
initially gapped state [l^, 32 1. For the latter the hori- 



zon effect is more pronounced as R is exponentially sup- 
pressed outside the light-cone. In contrast for our case, 
where the quench is from an initially gapless state, the 
change in the behavior of R from outside the light-cone 
to inside is less dramatic. 

Before we discuss the effects of the cosine potential 
Vsg on the correlator, we present results for the /3- 
function fil, m 



dg 
dlnl 



Kr, 



Kt. 



dK- 



d\nl 



dr] 
dlnl 



Ik{T, 



dlnl 



9 ? 

"■.9 7 

" 4 
ng'^j^K 
2 

2^,2 ] 



irg^^^K 



'^VTeff + -. I^ff (T„ 



dT„. 
dlnl 



(6) 

(7) 
(8) 
(9) 
(10) 



Eq. © shows that at long times after the quench (T™ ^ 
1), the periodic potential Vsg has a nonequilbrium scaling 



dimension given by Kneq- Since Keq < Kneq, the periodic 
potential is always more irrelevant for the nonequilibrium 
problem, however equations © show that the no- 
tion of irrelevance should not be taken literally out of 
equilibrium as the potential can cause inelastic scatter- 
ing that generates a dissipation rate 77, and at long-times 
a finite temperature Tej / Q . In this paper we will be in- 
terested in the pre-thermalized regime where the time is 
long as compared to microscopic time-scales, while short 
as compared to the dissipation-rate 1 <^ Tm ^ l/?7 
so that rj and Teff may be neglected. Since for weak 
quenches 77 ~ g^iKo — K)^ [§], this regime can be quite 
large. Eq. ^ also shows that even in the pre-thermalized 
regime, there is a crossover from an intermediate time 
dynamics where the physics is determined by the initial 
wave-function (and hence the initial Luttinger parameter 
Kq) and a long time time dynamics determined by Kneq- 
This can result in a situation where perturbation theory 



in g is violated at intermediate times when 



< 2 (i.e. 



g is a relevant perturbation in the initial state). We show 
below how this happens, and eventually present a RG 
treatment for the correlator in the regime where we are 
at weak-coupling at all times. In addition, we will study 
the system close to the critical point Kneq = 2+5, V(5 <C 1 
where logarithmic corrections to scaling arise. 

We evaluate R = +R'^^^ -f i?^^) perturbatively upto 
quadratic order in the periodic potential (i?^*^ being of 
0{g')). at a time after the quench is ^ 

dx' / dt' sgn{-c) 
Jq ^ 



(cos yj4>c{x't') - -4>{xiTra) - -4){x2T„ 



(11) 



The behavior of i?*-^' for several different quench proto- 
cols is shown in Fig. [T] We find that outside the light- 
cone, Tm < r/2, the correlator is non-zero. Moreover, for 

a quench where the initial state is such that the cosine 

2 TV- 
potential is relevant ( ^ ^ ° < 2), R can be parametri- 

cally large at these initial times outside the light-cone. 

When the initial state is one where the potential is irrel- 

2 Ty- 

evant or marginally irrelevant ( ^ ^ ° ^2), the correlator 
is most enhanced on the light-cone. For all these cases, 
within the light-cone Tm > r/2, the correlator reaches a 
steady-state. 

We derive the CS-like differential equation for R by 
splitting the fields 4> into slow and fast fields </) = </)< -|- <^> 
where the fast fields have a large weight at short wave- 
lengths, and therefore oscillate rapidly in time. We in- 
tegrate out the fast fields, and rescale the cut-off, po- 
sition and time. Such a procedure within the real-time 
Keldysh approach has been ernployed for the sine-Gordon 
model both for steady-state [8|, |9|, l33( and transient be- 
havior [l3|, where in each case the /? function was derived. 
Here we generalize this approach to the study of a two- 
point correlation function. This approach will reveal the 
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FIG. 1. The first order correction to the equal time correlation 
function i?'^'(r = 100, Tm) as a function of the time Tm after 
the quench. Length and time are measured in units of the 
cut-off A. Three different quenches are considered: a pure 
lattice quench {Kq — K = 2), lattice and interaction quench 
for {Ko = 1,K ^VS) and {Ko = 3,K ^VS). The light-cone 
is at Tm = r/2 = 50. We choose 7 = 2 and all quenches are 
to the critical point Kneq ~ 2. 



existence of a scaling regime where the correlator at long 
times and distances is related to the correlator at short 
times and distances (and a renormalized coupling), where 
the latter may be evaluated readily within perturbation 
theory. 

To leading order {g = 0) we find 



i?(°)(r,r„o-i?lc"^(r,r„) 



dA 

T 



(12) 



where i?< is the correlator of the slow fields, while R 
is the correlator for all the fields. The above implies 
the following differential equation (defining A = Ag/l, 
Aq being the bare cut-off, and Tmo is the time after 
the quench), - 7a„.o(r,r™)] ^kaAo) = q 

where (in units of A) 



7a«,o(?',T'„)= i 



K ' I 



2 il + (2r„+r)2 l + {2T„,-r)^ 



(13) 



Eq. (|13l) shows that there are two scaling limits where 
7an,o becomes independent of r,Tm- One is within the 
hght -cone Tm S> r ^ 1 where 7an,o — Kneq/2, and the 
second is on the light-cone r = 2Tm, r 3> 1 where 7aTi,o = 
{Kneq — Ktr/2) /2. Wc uow discuss the correction to R 
to next order in g where logarithmic corrections arise 
when Kneq 2 + (5 V(5 < 1. 

At next order, in terms of slow and fast fields, we 



find [29|, 

i?(i)(r,r™) =i?l3^(r,r™) 



+lanfi[r,Tm) — 



dA dA 

- 2 2K 

A A 

dA 



2g7rIcir,Tm)j^R%ir,t) (14) 



To quadratic order, the correction i?^^^ leads to the /3 
function which has been discussed above. Eqns. (|12p . 
(fT4| and the (3 function imply the following CS like dif- 
ferential equation for R, 



7an,o + 2TrgIc 



xR 



rAo AaTmo ' 







(15) 



where gi — g,S,r],Teff. We now discuss the solution to 
Eq. p5)) for two cases, one is for the pure lattice quench 
and the second for a simultaneous lattice and interaction 
quench. We will present results in the pre-thermalized 
regime Tmo <C l/?y where rj^Teff may be neglected. 

Pure lattice quench Kq = K and Keq = 2 + (5 V(5 <C 
1: Here we find. 



Icir,Tm) 



1 + 

1 + r2 



1 



- AT' 



4 - 8Ti 



4T2 



64r2 



(16) 



Note that Ic {Tm = 0) = as the lattice has not had time 
to affect the correlator. Eq. shows the appearance of 
scaling in two cases, one is within the light-cone where. 



(17) 



Ic{r, Tm » 1, Tm » r)= 1 + O ( -, — 



and the other is on the light-cone, 
Icir.Tm:^ l,2r„ = r) = 



(18) 



Eq. (|13p , ([TO)) also show that scaling is valid until (restor- 
ing units) I ^ Agr. In the scaling limit, Ik in Eq. (O is 
lK{Tm > 1) = TT 1-g^-f ... . Thus the solution 
of the CS equation dTSl) in the two scaling limits where 
ICi"fan,o, Ik are constants in time and position is. 



— rsC) -Ian (a') 

i?(Aor,AoT„,o,5o) = e J™ 
D I TmoAo , , 

x-R I — , — J — ,5(0 



(19) 



where 7a„ = 1 



— 2'nglc and the (3 function is 



dlfl = ^5*5' -Mri = -(27i"g)^. Note that the critical point 
corresponding to the = 1/2 Heisenberg chain corre- 
sponds to ^ = 2Trg. 
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The correlation function within the hght-cone is ob- 
tained by taking Tm ^ r, and integrating upto I = Apr, 



xR[ — = 1, — > 1,.9(0 



(20) 



where for ^ 3> 1, 5(0 being neghgibly small, the correla- 
tor on the r.h.s. may be evaluated within perturbation 
theory Moreover, = 1, 2^ > l) = 0(1) as this 

is a short-distance correlator. The solution of Eq. (^0)) 
at the critical point 5—2'iTg is found to be 



(Inr) 



/c-i/2 



r r 
In contrast the correlator on the light-cone is 



(21) 



i?(Aor,2r„,o = r,go) = e 
xi?('^ = l,2T™o = /,g(0 



(22) 



where the r.h.s is again a short distance correlator and 
therefore of 0(1). Due to Eq. (HH]), this leads to 



Rir 



Inr 

2^rn) — 

r 



(23) 



Eq. (pij) . ([25]) are the two main results of this section. 
Thus in the pre-thermalized regime the result Eq. (pij) 
within the light-cone for a pure lattice quench is the same 
as in the ground state of the Heisenberg chain [25-27|. In 
contrast, the quench leads to qualitatively new scaling be- 
havior for spatial separations on the light-cone (Eq. ([23]) ). 
Note that for a pure lattice quench, no dissipative effects 
are generated at long times to 0{g^) Q, extending the 
regime of validity of the pre-thermalized regime. 

Simultaneous lattice and interaction quench 
{Kq ^ K and K.^eq = 2 -I- (5,V(5 < 1): This quench 
corresponds to a final Hamiltonian whose ground state 
can be in the gapped phase. However due to the quench, 
since Kneq > Keq, the periodic potential is more irrele- 



vant (though it can give rise to inelastic scattering) [34 1 . 
Here we find that Ik in Eq. ([7]) is lK{Tm,Kneq = 
2)=7r [ci - C2 sin(7r(Xeq - Kneq)) lu (miu [Tm, l/rf\)] 
where ci,2 are 0(1) and depend on the initial Lut- 
tinger parameter Kq. When Kq—K, ci = 1. In 
what follows we will consider a prethermalized 
regime where T„i < l/rj. We will also assume, 
^1. In this case, Ik — tci is a 
The /? function in the vicinity of 



K, 



eq 



Kneq \ In T„ 



constant in time. 



K — 2 
where B = 



S becomes 



dg 
dlnl 

16 



X 3/4 



4 y::;7K^ - - In addition, Ic 
in Eq. (jl5|) becomes a universal function of the initial 
Luttinger parameter Kq and is plotted in Fig. [51 



iii 2- 




FIG. 2. Plot of Ic for a simultaneous lattice and interaction 
quench for points within the light-cone and near the critical 



point Kneq = 2 where Keq = ~' 



- 1. 



Using Eq. ([20)) to solve for the correlation function 
within the light cone, the correlation function is 



i?(Aor,r < T„o < l/v,9o) 



1 - (Aor)- 



1 



(Aor) 



l+A/2 



x/1 - (rAo)- 



2A 



1 + (Aor)-^ 



(24) 



where A = \/ — g^B'^ . Thus for an interaction and 
lattice quench, and for ^ Inr ^ 1, the correlators decay 
as a power-law with exponent (1 + A/2). This exponent 
is not the same as in a quadratic Luttinger liquid theory 
(1 -|- 5/2). In the vicinity of the nonequilibrium critical 
point A — )• 0, logarithmic corrections are obtained for 
A < Alnr < 1, where 



1 

i? - (Inr) « 
r 



(25) 



These corrections are significantly different from that 
near the equilibrium critical point 2^27|. In contrast 
to the regime within the light-cone discussed above, scal- 
ing behavior is lost on the light-cone for a simultaneous 
lattice and interaction quench because /c(r = 2Tm) ^ 



2 m 



r^t-rji .^yhere Ktr - 4 

In summary we have explicitly constructed a CS-like 
differential equation ^5)) for a two-point correlation func- 
tion for a quench in a spin-chain and found that for a 
lattice quench, scaling holds in two regimes. One is for 
spatial separations within the light-cone (Eq. ([HJ), and 
the second is for spatial separations on the light-cone 
(Eq. ([23| ). For a simultaneous lattice and interaction 
quench, the scaling on the light-cone is lost, but still holds 
within the light-cone (Eq. ([25]) ). The correlator is found 
to reach a quasi-steady-state where it decays in position 
with new power-laws that cannot be recovered within 
a purely Luttinger liquid model. These results hold in 
the pre-thermalized regime where dissipative effects are 
weak. For T„j > \/r\ dissipative effects will eventually 
make these correlators decay exponentially fast in posi- 
tion Q . How this crossover from power-law to exponen- 



5 



tial decay occurs in time, and observing it in numerical 
studies [35| are important open questions. 
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SUPPLEMENTARY MATERIAL: CORRELATION FUNCTIONS AFTER A QUANTUM QUENCH IN A 

SPIN-CHAIN IN THE PRE-THERMALIZED REGIME 

Model 



The model we will consider is that of an initial quench from Kq — >■ K, with the cosine-potential V = 
I dxcosljcj)) also switched on suddenly (at the same time as the quench). We use a bosonization prescription 
where, 



1/2 ^ 

L ^ \ 2t: J 



\\p\/2—ipx 



(26) 
(27) 



and A = a ^ is an ultra-violet cutoff. 
Thus the initial Hamiltonian for i < is, 

2ttJ 
= '^uo\p\vlvp 

while the time evolution from t > is due to 



■77- {dx(t>{x)f 



Hf = H 
Hfo 



V 



V 



f - ^/o 
u 

2^ 

gu 



dx 



K{TTn{x)Y 



{d^(j){x)f 



dx cos(70) 



(28) 



(29) 



(30) 



(31) 



Note that for L — >■ 00, the zero modes in Eqns. [26l [27]may be dropped. For finite-X, if one assumes that the quench 
connects the same zero-mode sectors of the initial and final Hamiltonian, they will again not play a role in the 
dynamics. 

To preserve Galilean invariance (which is not necessary for the formalism), we assume u = vf/K,uq = vf/Kq. 
Note that, 



cosh (3 — sinh (3 
- sinh (3 cosh /3 



7p 



cosh/3o -sinh^o\ f Vp 
-sinh/3o cosh/3o / [jj'Lp 



where e'"^^" = Kq, e"^^ = K. 
Let us define the functions 



f{pt)— cos(ii I p I t) cosh Po — i sm{u \ p \ t) cosh(2/3 — /3o) 
g{pt)— cos(u I p I t) sinh /3o + i sin(u | p | i) sinh(2/3 — (3q) 
which determine the time-evolution after the quench {t > 0) for g ~ 0, 

hl{t) + 6_p(i)- iript) - g{pt)) 7^1(0) + (fipt) ~ g*{pt)) 7y_p(0) 
hl{t) ^ b^p{t)= {r{pt)+g{pt))r^l{0) - U{pt)+9*{pt))il-p{Q) 

The Keldysh action is. 



Zk = Tr [p{t)] = Tr [ 
[<i)cu(pq\e 



-iHft 



e 

tiSo + Ssg) 



|V'j)('0»|e' 



IHftl 



(32) 
(33) 



(34) 
(35) 

(36) 
(37) 

(38) 
(39) 
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where So is the quadratic part which describes the physics in the absence of a periodic potential. In particular at a 
time t after the quench (note that the fields </> are real), 



'^0 = 7: [ dxi f dx2 [ dti ( dt2 (Ml) (t)q{i))(^-i 

^ J -00 J -co Jo Jo K'^R 

where 1 = (0:1,^2), 2 = {x2,t2) and 



G-/{1,2) 
(1,2) -[Gj,'GkG^'] (1,2) 



0c; (2) 



^cl,q = 



V2 



where 



and, 



Whereas, 



[G«-4(l,2)]-' = -S{x, - X2)5{t, - t2)^ [dl^,s - n'dl] 



0c; (1) 
0,(1) 



dx / dii[cos(70-(l)) -cos(70+(l))] 

-00 Jo 

For the quadratic theory after the quench, 

- mxt,)^{yt2))= r ^e-«f cos(p(x - y)) 
^ Jo P 

'u{ti -t2)+x-y 



Kq cos(u\p\ti) cos{u\p\t2) + -f^ sm{u\p\ti) sai{u\p\t2) 

^0 



K 



tan 



a 



+ tan"^ 



i{ti - t2) - {x-y) 



a 



The above implies. 



G^{xti,yt2) = -ie{h^t2){[cp{xti),<P{yt2)]) 
1 / u{ti - t2) + x-y 



= ~e{t,-t2) 



tan 



a 



+ tan 



1 / u{ti - t2) -{x-y) 



a 



G^{x -yM- t2) = G^'ix - y,-ih - ^2)) 
2 \ a 



tan 



1 / u{ti - t2) - {x-y) 



a 



and. 



G'^{xtt,yt2)= -i{{c^{xh),<j){yt2)}) 

Kq f dp _„„ . , ,, , I w 

= M+-^l J —e cos{up{x - y)) cos{u\p\{ti - t2)) 

Kq { dp , , , , 

-i— — e f cos(up(x-2/))cos(w|p|(ti +f2)) 



Let us define 



where a, 6 = ±. Therefore for 

C__(l,2) = e ^ V ' 
C++(l,2) = e M = 
C+_(l,2)=e" 





-t2)- 






-t2)- 






-t2)- 





(40) 



(41) 



(42) 



(43) 



(44) 



(45) 

(46) 
(47) 
(48) 

(49) 

(50) 

(51) 
(52) 

(53) 
(54) 

(55) 
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In general, 



g-^^[ "^-^2'"' + "^-1;'""' -'GA-(12)+mGA(l,2)+»bGfl(12)] 

It is convenient to define the following exponents that play an important role in the dynamics, 



7 



7 



Ko [1 

2 / ^2 



-^0 



Ko I 1 



^0 



In terms of these exponents, 

"2 ^iG^ {1,1) iG^ {2,2) 



11 
' 2 



iG'' {1,2) 



^/a^ + {u{ti-t2)+x-y}^ ^/a^ +{u{ti-t2) - {x ~ y)}^ 
in — h In — 



In, 



+ {u{ti + t2) + X - yy 
+ (2uti)2 



In, 



+ {u{h+t2)-{x-y)y 



(56) 

(57) 
(58) 
(59) 



(60) 



The above shows that the equal-time correlator for positions and times large as compared to the cut-off and for 
|r ± ut\/a ^ 1 is (setting a = 1), 



Gi{rt,Qt) = 







' 1 ' 


4 







(2ui)2 



r^2ut)^ 



-Ktr- 



(61) 



which agrees with Ref. |3lj who also point out that for 2ut <^ r, the correlator decays in position in the same way as 
in the initial state, but with a time-dependent prefactor (which is t^^*''). Exactly at the light-cone Ci{r = 2ut) ~ 

[M " [M "'■ ^o**^ th^* "i^i" t'^'^*' behavior of the correlator Rab = (e«'i'0a(i)/2g-»70t,(2)/2^^ which. 

follows in a straightforward way from that of Cab, has been discussed. 



Derivation of the /3 function from the action 



Let us denote a UV cutoff as 



A = 



We split the fields into slow {4>^) and fast {(jr') components [1, 0, [l^ 23 1 



6± = 0^ + 4>l 



and integrate out the fast fields perturbatively. In doing so to 0{g^), we obtain 
where 



(62) 

(63) 
(64) 



= / dR 



ut/\/2 ^ 

d{uT, 



2ttK 



Su 



(65) 



where Tg/ j — r\ — Q initially. 
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The correction from integrating the fast fields is, 



5S< = ^ r dx f dh [cos7(^<(l) -cos7<p<(l)] e'^ <K'''(^))'> 

U 7_co Jo 

pt pOO pt 2 

dxi / dh / dx2 / dt2e{ti - h) : cos (7<^<(1) - 7(/)<(2)) : e-^<(^-(i)-^-(2))^> 

Jo J—oo Jo 



oo 



l_e-^'<<^-(i)-^-(2)> 



(66) 



(67) 



e 2 



<(0+(l)-0+(2))^) 



1 _ e-7'{0$(i)0$(2)> 
2u2 



(68) 



/OO pZ poo pT 2 

dx, / dfi / da;2 / dt2{0{h-t2)+e{t2-h)} : cos (7</><(l) - 7</><(2)) : g" V<('^+(i)-'^-(2))'> 



where the fast correlators are related in a simple way to derivatives of the full correlators, 

dG 



G> =dA 



dA 



Explicit expressions for the fast correlators at equal time are 



2 



1 + 



-"■0 



2 



-"■0 



1 + {2tAy 



t»i/A dAKo 



1 + 



-"■0 



For the non-local fast correlators, we have 



(025(l)'^>/(2)) 



T 



4 



-"■0 



1 



1 



+ - 



^0 



l + A2((fi-t2) + (xi-X2)/«)2 l + A^{{ti-t2)-{xi-X2)M^ 

1 



l + Amh+t2) + {xi-X2)/u)^ l + Amh+t2)-{xi-X2)/uy 



6>(l),/.>(2)) = -z|^0(ii-f2) 



A((ii - fa) + {xi - X2)/u) 



+ 



A((ii - t2) - {Xi - X2)/u) 



(<^>(l)<A>(2)) = z|^^(i2-ti) 



l+A2((ti-t2) + (xi-,X2)Ai)2 l + A2((ti-f2)-(2;i-.T2)Ai)2 

A((ii - t2) + (xi - a;2)/u) , A((ti - ta) - (^i - .T2)/w) 



.l + A2((ii-t2) + (a;i-X2)/u)2 l + A2((fi-<2)-(a;i-X2)/u)2 
In the next step we define new variables 

_ Xi+X2 _ ti+t2 
r = Xi- X2 ,T = ti -t2 



Thus 



/CX) /-oo fOO p 

dxi / da;2 = / dR 
-oo J — oo J — oo J — 



dr 



dti / dt2 = 







oo J —oo 
dTjn I 

2T„ 



dT 



and perform a gradient expansion in R, Tm to obtain, 



5S< = 6S< + dS< + SS< + 6S< 



(69) 
(70) 

(71) 

(72) 

(73) 



(74) 
(75) 
(76) 

(77) 
(78) 

(79) 
(80) 

(81) 



where 



SS< = ^ r dx f dti [cos7(/)<(l) -cos70<(l)] e-^<('^c'W)'> 

J-oo Jo 



and 



OO f-OG 

dR / dr 

OO Jo 



2T™ 



-2T„ 



xlm 



-i-([0+(K+r/2,T„+T/2)-0_(ii-r/2,T„-r/2)]2) 



while 



^ig^A^dA 
^^T^ff - A 



OO /'OO 



4/72 



OO JO 



dR dr dTra dT{c^<{R,Tm)) 



2T„ 



2T„ 



xRe 



e 2 



:L-([0+(iJ+r/2,T^+r/2)-<^_(fl-r/2,T^-r/2)] = 



and 



.c<_ gW c;A 



/OO /.OO rt/V2 i'2Tm 

dR dr dT^ / dT<P<{R,Tm)TdT^<P<i{R,T, 

-OO J-oo ^0 J -2Tm 



xlm 



where 



1 



1 + A2 (r + r/u)^ 1 + A2 (r - r/uf 



tr 



1 + A2 (2T™ + r/w)' 1 + A2 (2T™ - r/u)' 



-iK, 



eg 



A(t + r/u) 



A(t — r/ii) 



1 + A2 (r + r/w)^ 1 + A2 (r - r/uf 



whereas 



e 2 



{[4>+{R+r/2,Tr„.+T/2)-,l,_{R-r/2,Tr„.-T/2)f) ^ 



v/1 + A2(t + r/M)2 ^1 + A2(t - r/w)2 



y/l + A2{2(r„, + t/2)}2 ^1 + A2{2(r^ - t/2)P 
Vl + A2(2T„, + r/«)2 ^l + A2(2T„-r/M)2 

^g-i/fe<,[tan-i(A(T+r/«))+tan-i(A(T-r/M))] 

where i?e[A] = (A + A*)/2, /m[y4] = {A-A*)/{2i). 
Collecting all terms we find, 



d{R/u) / dT^ [cosjcP<{R,Trr,) - COSJ(t><{R,Trr,)] e-T(K=l(T^)) > 

-OO Jo 

/OO pt/ y/2 

d{R/u) / dTm [-lR{Tm) {dR/u4>tl) {dR/u<t't) " (^^-) (St^<^c1) (5t^<^9 )] 

-OO JO 



2„,2 J A /-OO 



55n< = 



^2^2 



° " 2 A 
< _ig^^f^dK 



5S< 



gfj^AdA 
2~T 



rt/V2 

Jo 
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where 



MTm) = -A* / d{r/u) / dT{r/uflm 



2T„ 



d{r/u) / rfrr^/m 

-oo ^0 



2T„ 



2T™ 
2T„ 



drRe 



d{r/u) I dTTim 

-oo J-2Tjn 



^-\{[4>+(R+T/1,T^+Tl1)-,t>-(R-T/2,T^-T/2)Y)pl^^^rp^^^ 

^-^{[4>+{R+r/2,T^+T/2)-4>-{R-r/2,T^-T/2)f)p^^^j,^^^'^ 
^-^{l<l>+{R+r/2,T„+T/2)-<l>_(R-r/2,T^-T/2)f)p^^^ j,^^ 
^-^{[4>^(R+rl2,Tm+r/2)-<t>-(R-rl2,Tm-T/2)\^)p(^^^j,^^^' 



Ik ~ Ir — It^ 
Iu= Ir + 



(92) 
(93) 

(94) 

(95) 

(96) 
(97) 



At the next step we rescale the cut-off back to the original value of A, and in the process rescale position and time 
to R,Tjn ^{R,Tm)- This rescaling is not necessary in expressions for SS^j.^^^^^ as they are already of O (^)- We 
also express everything in dimensionless units of ^ = AR/u,T = TA. Thus to summarize one obtains, 



: / dR [ 

J -oo Jo 



-^2 , q2 \ . < 



and 



A 



6S< 



A' 

gW dA 



dR 

oo ^0 



^^T^tf - 2 A 



5S< 



g^dA 
2 A 



dfm [cos70<(i?,T„) - cos7</'^(i?,T„)] e-^<(^Ji(^'"))'> 

/oo rth/yj^ 
^ dR df^ [-lR{Tm) {dn<t>ti) {dR<t>t) - It^ (Tm) {df^<i>fi) {df^<t>t)] 

/ dR / df^{4't) lT.,,{Tr.) 

J-oo Jo 

tA/V2 



dR 

oo Jo 



Since, 



, dA , A - A' 



dln(/) 



' A ' A 

Therefore the RG equations in terms of dimensionless variables such as Tm — > TmA, ry — ^ jj/A, Tg// — ^ Te///A are 



(98) 

(99) 
(100) 
(101) 
(102) 

(103) 

(104) 



dg 

d\nl 



9 



"^9 ' 1 + 47^2 

IxiTm) 



dK-'^ _ ng^"f^ 
d\nl ~ 4 
1 du 7r(/^7^ 



dlnl 



Kudlnl 4 
df] ng'^j^K 
rfh^ = ^+^^'''(^'"^ 
ng'^^'^K 



= '^VTeff + 



4 '^'ff 
dTr 



dlnl 



T 



(105) 

(106) 
(107) 
(108) 
(109) 
(110) 
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Note that the renormahzation of the velocity in Eq. 11071 is a minor effect which wih be neglected in the main text. 
To summarize: 



/fl(T™A) = - 
/(,(T„A) = - 
/t,„(T„A) = 
/„(T„A) = - 
/„(T„,A) = 



dfr^Im [B(f,T„A,f)] 
dff^Im [B{f,T^A,f)] 



-co Jo 

p2T,„A 
df / ^^^^"^1 

-oo Jo 

/2T„A 
dfRe[B{f,T^A,f)] 
-2T,„A 
oo /•2T„A 

df / dfflm[B{f,T,nA,f)] 



2T„A 



oo 

df 

-oo Jo 



2T„A 



df(f2 +f2)lm [B(f,T,„A,f)] 



/oo /•2T„A 
/ df (f2 - f2) Im [B(f , T„A, f)] 
-oo Jo 



lK{TmA) 



where Re[B] = {B + B*)/2, Im[B]={B - B*)/ {2i), and 

B(f , T„A, f) - C+_ (f , T„A, f)F(f , T„A, f) 



(111) 
(112) 
(113) 
(114) 
(115) 
(116) 

(117) 



with C+_(f,T„A,f) = (^pt'P+{f,T+T^A/2) ^-ici>^{o,f-T^A/2)y rpj^jg quantity within leading order in perturbation theory 
is, 



C+_(f,T„A,f) 



1 



1 



v/l + (f + f)2 Vl + (f-f)2 



while F is given by 



y/l + {2(T,„A + f/2)}2 ^1 + {2(r„A - f/2)p 
v/1 + (2T„A + f)2 ^1 + (2T„A - f)2 

^ g-iife, (T+f )+tan-i (r-f )] 



F{r,TmA,T) = Kneq 
1 



-ft:* 



(118) 



1 + (f + f)^ 1 + (f - ff 



1 + (2r„A + ff 1 + (2r„A - ff 



T A- r 



1 + (f + f)^ 1 + (f - ff 



(119) 



Perturbative evaluation of equal time correlation function 

Let us study the following two-point function to next to leading order in perturbation theory, where 

i?(x,^,x.^)^4(cos(^)cos(2:^)) 



where 



f(12)] 



(120) 
(121) 

(122) 
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Using Eq. |60] we get 



2K„ 



1 



1 



1/4 



(123) 



The above may be rewritten to give in the scahng hmit, Eq. 5 in the main text. 
To next order the equal time correlator is, 



R'-^Hxit,X2t) =^-^ [ dx' f di' V s5n(-c)4(cos(7(?!)c(a;'t'))cos (^0(a;it)) cos (^0(a;2i))) (124) 



c=± 



On evaluating the expectation value, one finds, 

'-/"^K ( _^ (u{t-t')+xi-x' 



dx' I dt' sin 

-OO Jo 

-t'^K ( _i / u{t - t') + X2-x' 



- <|^tan 



tan 



1 f u{t- t') - {xi - x') 



a 



-■j^tan 



a 



tan 



1 ( U{t - t') - {X2 ~ x') 



(125) 



The above expression for i?'^^^ for different quench protocols is discussed and plotted in Fig. 1 of the main text. 
At short times 



i?W(xi,..;tA«l)^g('2:^) {Mf 



1 



1 



i?(l)(xi,X2;tA« l)^^l^^g ( ff) (At)^i?(")(|xi -X2|)ci 



, \xi - a;2| > 1 



(126) 
(127) 



where a [6, \xi - X2I] = dx j^^^2]i+t.[i+(^_(^^_^^))2]« is independent of - X2\ for - xs] > 1 



Derivation of a Callan-Symanzik-like differential equation out of equilibrium 

In order to derive a CS-like differential equation for the equal-time correlation function 

Rix,t,X2t)^4{cosi:^] cos^^^("^') 
= + + . . . 



(128) 
(129) 



we split the fields (j), into slow and fast fields, and integrate over the fast fields. In doing so, for example, i?^"^ may be 
written as, 



K„ 



(xi-:c2)"A"/»" Ktr Ktr 

+ (xi-ic2)^A^/u^ l + (2tA)^ 2 



R'^°\xit,X2t) ^ R%{xit,X2t)e 

n we rescale the cut-off, and 1 
Expanding the above expression to O (^), and noting that 

A 



71+ 



l + A'i(2t+(cci-x2)/ti)^ l + A'^(2t-(xi-x2)/i»)^ 



(130) 



Then we rescale the cut-off, and also rescale position and time. In doing so, the above expression does not change. 



(A) - i?!;") (A - dK) ^ 



dk 



A- 



(131) 
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we obtain the CS-like equation to leading order, 

d 



(9 In A 

d 

dhTl 



where Aq is the bare cut-off and, 



fi{xit,X2t) 

(Xi - X2)2AV«2 



' 2 U 



l + {xi- X2Yh?/u^ 1 + (2tA)2 

1 



+ 



+ A?{2t + (a;i - X2)/uY 1 + A2(2i - - X2)/u 



Some limiting expressions for 7an,o are 



1 'y'^Kn 

7an,o(|a;i -X2|A> l,tA< 1) = -{Kneq+Ktr) = "^-^ 



where 



S{xi,X2,t,x' ,t') = 



2i^„e, + ^ ^ A)2 ' l + (2iA)2 



+ 



2 [l + A'^{t-t' + {xi-x')/uf l + A2(f-t 

Ktril 1 

2 



- {X, - X')/Uf } 



l + A2(f + t' + (.Ti -,t;')/w)2 1 + A2(t + /:' - (xi -.t')/m)2 

j^ne, [ 1 1 

2 \l+h?{t-t' + {X2-X')/Uf l+h?{t-t' -{X2-X')/UY 

Ktr ( 1 1 



2 [l + A2(t + t' + (a;2-a;')/u)2 1 + A2(t + - (a;2 - a;')/u)2 
.^Keq( K{t-t' + {xi-x')/u) K{t-t! -{xi-x')/u) 

2~\l + A2(i-f' + {xi - x')/uY ^ l + K'^{t-f' - {xi-x')/uY 

Keq r A(i ~t'+ {X2 - x')/u) 



-ic 



-ic- 



2~\l + A2(i-i' + ( 



A{t-t' - {x2 -x')/u) 

X2-X')luf l+K^{t-t' -{X2-X')/Uy 



(132) 



{\xi — X2\A » l,tA » l,2t » |a;i — X2\) = 
7an,o(|a;i - a;2|A > l,tA > l,2t «: \xi - X2\) = 

1 / K 

7an,o(|a;i - a;2|A > l,iA > l,2i = \xi - X2\) = - ( if„eg ^ 

Now for we find, 

Ri'\x,t,X2t) = 4TV p(^o<+^o>)i5.,(<^< +0>)cos (^7^!(M±2^1(M) eos (^ 7^<M) + 7.^>(x20 y 
After some algebra we obtain, 

R^^\xit,X2t)=^-^ f dx' I dt'^sgn{-c) 

^ J —oo Jo c=± 

(cos (70<(^'i') - I'l^^ixit) - ^</'<(a^2i))) 
x(cos {jcl>>{x't') - ^e{x^t) - ^-e{x2t))) 

J — oo Jo c=± 

x(cos (7<^<(xY) - |<^<(a;if) - |<^<(a;2i))) 



(133) 

(134) 
(135) 
(136) 
(137) 

(138) 



(139) 



(140) 



(141) 
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Rescaling the cut-ofF and correspondingly the position and time, we get 



. ^dA dA 
1 + 2— - 2K, 



dA 



{Xi - X2,t) — 



dA r, 

-2gnIc(xi,X2,t) — R^{xi - X2,t) 



Rewriting the zero-order term for convenience, 

R^°^{xi~X2,t)=R'^^\xi-X2,t) 

and combining Eq. 11421 11431 we get 



dA , ' 

1 - -J-^an.a(^l - X2,t) 



(142) 



(143) 



R^ R 



(0) 



dA ^ ^ dA^ ^ , , 

1 —lanfi\X\ - X2,t) + —2g'Klc[Xi - X2,t) 



A 



. dA 



I ^ I dA 



(144) 



The above imphes. 



d 



+ lanfi - 2TTgIc 



R{xi,X2,Tm) = 



91nA 

Including the /3- function derived in a previous section and which appears at 0{g'^), the above can be rewritten as 

VAo AqT^o 



(145) 



R 



= 



(146) 



where gi = g, K, 77, T^ff, and Ao, T„ia,gio denote bare values. Moreover Ic is given by 
Icixi,X2,t,K^,,,K,,) ^ ~ (J-^ e^"-'"Vi+(-i-^)^ dx' J^dt' 



sm 



xe 



K ( 

j tan"^ {t' + x' - xi) + tan"^ (t' - {x' - xi)) + tan"^ (t' + x' - X2) + tan"^ [t' - {x' - X2)) 

[in ^l + (t'+a:'-i;i)='+ln ^l + (t'-(x'-xi))2+ln l + (t' +x' -Xif +\n ^l + (t'-(i:'-a:2))^] 



In 



l + (2t-t' -iciY 



+ln 



- + (2t-t' -ii)Y 



l+4(t-t')^ 



+ln 



hln 



l + [2t-t' -{31' -^2)) 
l + 4(f-t')-^ 



neq 



1 + 4^2 ' l + 4(i-i')2 2 ll + (t' + a;'-a;i)^ 1 + (<' - (a;' - xi))^ 1 + (<' + a;' - a;2)^ 



1 



1 



1 



1 



1 + ~ {x' ~ X2)f } 2 {l + {2t-t' + x' -xi)^ 1 + {2t - t' - {x' ~ xi))^ l + {2t-t' + x' -X2y 
1 



1 + (2t - i' - (X' - .X2))2 

+— e^"-'"Vi+(-i-2)^ dx' / dt'cos [^|tan-i(i' + a;'-a;i)+tan-i(i'- (x'-xi)) 
27r J_oo Jo L 2 L 

+ tan"^ {t' + x' - X2)+ tan"^ (t' - {x' - X2)) 



(147) 



xe 



[in y/l + {t'+x'-xiy^+ln ^l + (t'-(x'-xi))2+ln Vl + (t'+a;'-a;2)^+ln Vl + (t'-(x'-X2))=] 



In 



l + (2t-t'+x'-xi)^ 



+ln 



if. 



l+(2f-f'-(x'-xi))" 
l+4(t-t')^ 

(f - (a;' - xi)) 



+ln 



l+(2t-t'+a;'-X2) 



hln 



l + (2t-t'-(x'-j:2))'^ 
l + 4(t-t')^ 

(<' - (a;' - X2)) 



(t' + x'-Xi) ^ (t'-jx'-Xi)) ^ {t' + x'-X2) 

l+{t'+x'-Xif 1 + {f - {x' - Xi)f l + (t' +X' -X2f^ l + {t' -{X' -X2)f 



(148) 
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The usual CS equation encountered in equilibrium for the Heisenberg chain 2^27| may be obtained from the above 
by setting the time after the quench t = oo, and by setting Kq = K and by also in addition being at the critical point 



Keq — 2. Here, noting that — Im 



g — ztan ^ X 1 

l + ix 



l+x'-' 



1 + 



^ cos (tan ^ x) , one may write 



(1 + (X1-X2)2) 

Ic = Im 



2n 



oo /"OO 

dx' / (it'e~ ln(l+^(t'+^'-^i))-ln(l+^(*'-^'+^i))-ln(l+^(t'+^'-^2))-ln(l+i(t'-a:'+a;2)) 



1 



-oo JO 
1 



+ 



1 



{l + i{t' + x' - xi) 1 + i{t' - x' + xi) I + i{t' + x' ~ X2) I + i{t' - x' + X2) 
Analytically continuing in time, 



00 Jo 



d_ 



(l + r)2 + (x'-xi)2 (l + r)2 + (x'-X2)2 



^ (xi - x^f + 1 
4+ - X2)2 

Thus for |xi - x-2\ > 1, 



(149) 

(150) 
(151) 

(152) 



Correlation function for the pure lattice quench 



We want to solve 



^+/3(5.)|--7a„,o + 2.g/c 



R 



(xi - a;2)Ao Ao^, 



mO 



I 



(153) 



where gi = g, K,r],Teff. In addition, Ao,Tmo, gio denote bare values. In this section, we want to solve Eg. 11531 for 
the lattice quench Kq — K and near the critical point K^q = 2 + 5. Note that. 



7a«,o(|a;i -2:2|A> l,tA< 1) = - [Kneq+Ktr) = 



l^Ko _ S 
8 2 



7a«,o(|a;i -X2|A> 1,<A> l,2t > -a;2|) 
7a«,o(|a;i -a:2|A> l,tA> l,2t < -a;2|) 
7a«,o(|a;i - X2IA > l,tA > l,2t = |a;i - X2I 
and for Ic , the full expression in Eqns 11471 11481 reduces to 



^neq ^ 

2 ^ 2 

2 ^ 2 
--i K - — 



= 1 



(154) 
(155) 
(156) 
(157) 



(Xi-X2)^ + 1 (^1-^2)^, 
ic(a;i - X2,t)= -, z Im 



A+{xi-X2y 



2tt 



00 

I I dx 



1 



1 



(l + ii)2+a;2 (I + ^tY + {^j: - {xi - X2)f 



(Xi -X2)2 + 1 _ 1 + (.Ti - X2f 

4+ (xi - X2Y 



l + t2 



{xi - X2f - 4i2 + 4 - 8i2 
(a;i-a;2)2-4t2 + 4l + 64i2 



(158) 



The following different limits exist, 

Ic{\xi - X2I > l,t > l,t > \xi - X2\)= l + O 



[xi - 2:2)2 ■ 



{Xi - X2f 
t4 



Ic{\xi -X2I > l,t > l,t< |a:i -2:21)= 1 + ( ^ 



/c(|a;i - a;2| > l,t > 1, 2t = [xi - X2\)= 1 + 1/2 = 3/2 + O 



1 



{xi - X2y 



(159) 
(160) 
(161) 
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In addition the /3-function for a pure lattice quench and using Kn^q = Kf,q = 2 + (5 is, 

dg 



dS _ 
dhd ~ 
1 du 



udlnl 
dr) 



dlnl 



9^ 



J] + Ang IniTrn) 



d\nl 

dT„, 

1 rr 



d\nl 



where 



thus 



lK{Tm, Kneq = Keq = 2) = TT — — 



2 + 20Tl + 7 X 16T^ 

(i + 4r2)3 



lK{Trn < 1) - 2ttT^ 

7 



lK{Tm:> 1) = 7r 



1 - 



Moreover, 



(2r,„)' 



Ir^ff = 67rT„ 



[1 + (2T„02 
1 - ^T^ 

3 rn 



.(1 + 4^2)3 

At very long times, for the pure lattice quench, we may neglect Teff,r] as these are zero at steady-state at 0{g^) 



(162) 
(163) 
(164) 
(165) 
(166) 
(167) 

(168) 

(169) 
(170) 

(171) 
(172) 



Now, Eq. 11531 implies 



where 



Eq. 11731 implies 



R 



rAo T^qAo 



I 



I 



g{l) ] =e^«o"'^s'^^i?(Aor,r,„oAo,ffo) 



Ian = 7an,0 " 27rg/c 



i?(Aor,AoT,„o,<?o) =e ^™ "'^ WFi? 



I ' I 



,3(0 



(173) 



(174) 



(175) 



One interesting case is the behavior of correlators at steady-state Tmo oo. Here the explicit dependence on the 
time after the quench drops out. Then integrating upto / = Agr, we get 



i?(Aor,r,„o = oo,.9o) = e ''a i? f _il = l, T„,o - (X3, g(0 



(176) 



where for I ^ 1, g(X) being negligibly small, the correlator on the r.h.s. may be evaluated within perturbation theory. 
Moreover, R (^^ = l,rmO — oo,g{l)) = 0(1) + O {g{l)) + ... as this is a short-distance correlator. 
Another interesting case is for spatial separations on the light-cone where r = 2Tmo. Here, 



R (Aor, 2T„iO = f , .go) = e 



dg'- 



xR[^ = l,2T„,o^ 1,9(1) 



(177) 
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Here R on the r.h.s may again be evaluated within perturbation theory, and being a short-distance correlator, R 
Oil). 

At the critical point S — 2-Kg, so that 



d5 
dhi7 



(178) 



and 



7an = 1 + 2 ^ ^^C 



(179) 



where /p = 1 for T,„o S> whereas Ic — 3/2 for points on the light cone r = 2TrnQ- Thus in general. 



1 1 



13 (52 g 



(180) 



The solution of Eq. [TZ8]is 



^ - ^ = In (Aor) 



(5(Aor) 6o 



(181) 



while 



fSiAor) / 1\ 

1 .i^f ^...M + (fc-5).- 



So 



(182) 



Thus the correlation function is found to be 



R (Aor, go) 



[In(Aor)] 



/c-l/2 



Anr 



(183) 



Eq. 11831 shows that far from the light cones where T^o ^ r, the correlator decays as as in equilibrium, while 

right at the light-cone r = 2Tmo, the correlator decays somewhat slower, as — . 



Correlation function for the lattice and interaction quench 

Here we are interested in the regime where K„i,q ~ 2, which implies 



Ktr - - 2 (185) 

Writing Kneq = 2 + 5, the RG equations for times T,„ < 1/?/ (i.e., on neglecting dissipative effects) are the following 
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above we have used that since Kneq = 2 + S, dK = ^^dS, and K^q = ^ ^° a /^j^ ^ 1 ■ Moreover, 



/>oo 

IxiTm = OO, Kneq = 2) = -4 / du 

Jo 



du COS (Ken tan ^ m) 

1 + 



(1+U2) 



dv- 



COS (iiTeq tan m) 



dv- 



1 + v 



■ sin (i^eg tan^i v) (189) 



(l+t;2)2 



sin {Keq tan ^ w) 



-2K 



eg 



i2(M) 
2 



f du r- cos ( /Ten tan u) / dv- r-^ cos ( tan v) 

1 + ' Jo (1 + ^ ) 



-2Ke 



du 



(1+^^2)2 



fU 

sin fi^eo tan^^ u) / :r sin (Ken tan^^ ' 



(190) 



(191) 



(192) 



Note that Eqns 11901 11911 are logarithmically divergent (such a divergence does not arise for a pure lattice quench) . 
The reason for this divergence is because we have used leading order perturbation theory to evaluate the correlators 
that go into Ik, whereas, we should have used the correlators from the full theory which correspond to a non-zero 
dissipation 77. Taking this into account, the logarthmic divergence is cut-off by min(rm, I/??)- Denoting the well 

behaved terms in Eqns fT89l 11921 as ttci, and using 12(00) = — -4, ^in { ^^2"' ) ^^"-^ *3(oo) — k^{k'^ -4) sin ^ ^^^'' ^ , 
we may write, 



IxiTm^Kneq = 2) = TT 



ci — 2c2 sin 



ttK, 



eq 



min{T^,l/rj) 



du- 



U 



l + U 



2 COS (^Keq tan ^ u) 



TT [ci - C2 sin (Tri^eq) lu (min(T„i, 1/??))] 



(193) 
(194) 



where ci^2 are 0(1) and depend on Keq and hence Kq. When Kq — K, ci = 1. 

In what follows we will consider a prethermalized regime where Tm < l/rj, and \Keq — Kneql^nTm ^1. In this 
case, Ik — ttci. The solution of the RG equations in this regime, under conditions that 



9eff ^gB <S 



where 



B = -K^/ci 



j'^Ko\ f 16 



72/10 



3/4 



(195) 



(196) 



and defining a constant of the RG flow -jf^ — —geffS, -^^^ = ~giff 



A^JP-g, 



is given by. 



5{l) = Acoth 
9eff{l) = 



A\nl + tanh" 
A 



-1 M 
.So 



sinh 



Aln? + tanh"^ ( # 



(197) 

(198) 
(199) 



We will now study the system in the vicinity of the non-equilibrium critical point Kneq — 2 + 6. Here for Iq we find 

2' 



/c(r> 1,T,„> 1,T„. >r 



Ke. 



1 



Ke. 



L(Keq) 



(200) 
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where 



L{K,g) = lim f 
|r|»iV47r7o 



2 pOO r-U 

du I dv sin 



^-3 J 



tan u + tan d + tan {u + ?-) + tan (v — ?') 



VI + m2 Vl + ^1 + (u + r)2 ^1 + - r)2 

+ (r ^ -r) 



1 + 1*2 1 + ^,2 1 + (u + r)2 1 + (u - r)2 



(201) 



L(Keq) is a universal function of iiTeg and hence the initial Luttinger parameter. Ic is plotted in the main text. 
The solution of the CZ equations is, 

i?(Aor,Aor„,.o,5o) = e ^ ^hTTRI—, — ,g{l) 



(202) 



We will solve this for l/rj ^ Tmo ^ r. For this case, dissipative effects may be neglected in the solution of the RG 
equations which will be integrated upto / ~ Aor. In particular. 



R 



(Aor,r « T,„o « l/v,9o) = e"^™""" ^^'^i? = 1> ^ » l>5(o) 



(203) 



In the r.h.s, R being a short-distance correlator, is 0(1). Note that 



(204) 



Then, 



9(0 ^(a') 1 



sinh ( y4 In ^ + tanh' 



-1 A_ 
So 



sinh ( tanh ^ ^ 



27r/cA 
B 



In 



tanhf 4ln; + itanh"^ # 



5o 



tanh tanh ^ 



(205) 



3/4 



wherei3 = ^V^(2^) - 1^ 

The correlation function within the light-cone is given by setting I — r. 



i?(Aor,r < T„o < l/fy,5o) 



(Aor)i+^/^V 

For Alnr ^ 1, we obtain the following logarithmic correction to scaling. 



1 - (Aor)- 



1 + (Aor)-^ 



1 2irle 1 

i? - (Inr) « 2 
r 



(206) 



(207) 



Let us now discuss how Ic behaves for spatial separations on the light-cone. Simplifying the full expression in 
Eqns. [TiTlfTigl bv noting that for r > 1, t > 1, r = 2i, 



'i + {2t-t' + x'y 



1 + 4<2 



'l + (2i-<'-a;') 
l + 4(<-t')2 



\ + {2t-t' + x' -rf l + {t'^x'y 



1 + 4t2 



1 + 4i2 



'l + {2t-t' -x' +rY 

A2 



1 + 4(t - t') 
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one may write, 



/c(r = 2t, r > 1, = 2) = 2-^*-/2Lt, 



27r 



dx' / dt' 

-OO Jo 



Sin 



1 



tan"^ {i! + x') + tan"^ [t' - x') + tan"^ [t' +x' -r)+ tan" 
1 



-i(i'-a;' + r)| 



^Jl + {t' + x' -rf ^Jl + {t' -X' +rf 



1 



(1 + ^) 



l + (t' + a;')^ l + {t'~x'f l + {t' + x'-rf l + {t'-x'+rf 



'K. 



2-K 



r/2 

dx' I dt' 



OO Jo 



COS 



jtan"^ {t' + x') + tan"^ {t' - x') + tan"^ {t' + x' - r) + tan"^ {t' - a;' + r)| 



1+ 



Y^l + (t' + x')^ VV'l + li'-a;')^ 
(i' + x') (f - x') 



1 + + x' -rf Jl + [f ~x' + rf 



[f + x' - r) {f -x'+r)) 



+ 



l + {t' + x'f l + {t'-x'f l + {t' + x'-rf ' l + {t'-x' + rf 



Here we find that 



Ic{r > 1) r 



(208) 



(209) 



Thus the scaHng for spatial separations on the light-cone r = 2t which we found for a pure lattice quench, is lost for 
a simultaneous lattice and interaction quench when Ktr > (i.e., > 2 so that the initial state is irrelevant). In 
contrast, scaling continues to hold for spatial separations inside the light-cone. 

I 



